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Engineering flat electronic bands in quasiperiodic and fractal loop geometries 
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Exact construction of one electron eigenstates with flat, non-dispersive bands, and localized over 
clusters of various sizes is reported for a class of quasi-one dimensional looped networks. Quasiperi¬ 
odic Fibonacci and Berber fractal geometries are embedded in the arms of the loop threaded by a 
uniform magnetic flux. We work out an analytical scheme to unravel the localized single particle 
states pinned at various atomic sites or over clusters of them. The magnetic field is varied to control, 
in a subtle way, the extent of localization and the location of the flat band states in energy space. In 
addition to this we show that, an appropriate tuning of the field can lead to a re-entrant behavior 
of the effective mass of the electron in a band, with a periodic flip in its sign. 
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Geometrically frustrated lattices (GFL) supporting 
flat, dispersionless bands in their energy spectrum with 
macroscopically degenerate eigenstates have drawn great 
iirterest in recent tiirres [l|-@- Initial interest in antiferro¬ 
magnetic Heisenberg model on frustrated lattices [^-fl^ 
has evolved into extensive studies of the gapped flat band 
states to gapless chiral modes in graphenes [l^, in opti¬ 
cal lattices of ultracold atoms [l^ , waveguide a rray s [l7l| , 
or in microcavities having exciton-polaritons |18l| . The 
quenched kinetic energy of an electron in a flat band state 
(FBS) leads to the possibility of achieving strongly corre¬ 
lated electronic states, topologically ordered phases, such 
as the lattice versions of fractional quantum Hall states 
M- Recently, the controlled growth of artificial lattices 
with complications such as in the kagome class has added 
excitement to such studies [l^- plj . 

Spinless fermions are easily trapped in flat bands 
The non-dispersive character of the energy (E)- wave 
vector (fe) curve implies an infinite effective mass of the 
electron, leading to practically its immobility in the lat¬ 
tice. Such states are therefore strictly localized either on 
special sets of vertices, or in a finite cluster of atomic 
sites spanning finite areas of the underlying lattice. Re¬ 
cently it has been shown that, an infinity of such cluster- 
localized single particle states can be exactly constructed 
even in a class of deterministic fractals [l^. Apart from 
its interest in direct relation to the study of GFL’s, this 
work provides an example where eigenvalues correspond¬ 
ing to localized eigenstates in an infinite fractal geometry 
can be exactly evaluated, a task that is a non-trivial one 
if one remembers that these fractal systems are free from 
translational invariance of any kind. 

In this communication we unravel and analyze groups 
of flat, dispersionless energy bands in some tailor made 
GFL’s. The lattices display an interesting competition 
between long range translational order along the horizon¬ 
tal {x—) axis and an aperiodic growth in the transverse 
directions. In each case, the skeleton is an infinite array 
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FIG. 1: (Color online), (a) A typical realization of an elemen¬ 
tary array of diamond shaped quadruplets exhibiting a single 
flat, dispersionless band at A = 0, (b) its renormalization into 
an effective two-arm ladder network with energy dependent 
potential and nearest and next nearest neighbor interactions 
and, (c) the effective one dimensional chain obtained after 
decimating out the top (yellow colored) vertices. 


of diamond shaped networks (Fig. [TJa)) where hierar¬ 
chical structures, grown deterministically, are embedded 
in each arm of the diamond. A uniform magnetic flux 
threads each elementary plaquette, as will be illustrated 
in appropriate cases. The motivation behind the present 
study is two fold. 

First, exploiting the self similar growth of the hier¬ 
archical structures we can implement a real space renor¬ 
malization group (RSRG) scheme to evaluate the pinned, 
flat band states (FBS) exactly. In all the cases we discuss 
in this paper the flat band states merge with the local¬ 
ized eigenstates of the infinite hierarchical geometry as 
higher and higher generations of them are embedded in 
the arms of the underlying diamonds. Thus, once again, 
the problem of an exact evaluation of localized states in 
a fractal geometry is, at least partially, solved. 

Second, and a more important aspect of the problem 
is to look for a possible coexistence of a multitude of 
dispersive and non-dispersive energy bands in a periodic 
array of the elementary building blocks. As an arm of 
the elementary diamond hosts a quasi-periodic or a de- 
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terministic fractal of sequentially increasing generation, 
the flat, dispersive bands are likely get densely packed 
in an environment of dispersive ones. The density of 
packing may even lead, for a large enough generation of 
the hierarchical structure, to a re-entrant dispersive to 
non-dispersive crossover in the behavior of the electrons. 
A recent inspiring work by Danieli et al. @ has shown 
that a quasiperiodically modulated flat band geometry 
may even allow for a precise engineering of the mobility 
edges. 

In addition to this, we anticipate that, a variation in 
the strength of the magnetic field can lead to a tun¬ 
ing of the curvature of the energy dispersion curves. 
We can thus achieve a comprehensive control over the 
group velocity and effective mass of the electron with the 
help of an external field. Knowing that, a deterministic 
quasiperiodic or fractal geometry normally offers a com¬ 
pletely fragmented. Cantor set-like energy spectrum, this 
latter study may allow us to control the effective mass of 
the electron using an external agent such as the mag¬ 
netic field over arbitrarily small scales of the energy or, 
equivalently, the wave vector. 

We present here a simple analytical method to detect 
the sharply localized eigenstates that are pinned on cer¬ 
tain atoms or atomic clusters in a periodic array of di¬ 
amonds. The non-dispersive character of such states is 
explicitly worked out. The method is then extended to 
unravel the entire set of such FBS when each arm of an 
elementary diamond hosts a quasiperiodic latti ce g rown 
according to a deterministic Fibonacci sequence (23|, and 
Berker |24| geometry. Such aperiodic structures with se¬ 
quentially increasing hierarchy are embedded in the dia¬ 
mond’s arm. We indeed find that, as one gradually in¬ 
creases the degree of aperiodicity in the arms, the pinned 
FBS in such periodic approximants turn out to be the 
localized eigenstates of the systems in their respective 
thermodynamic limits. 

In addition, we observe that, with a deterministic ape¬ 
riodic geometry of sufficiently large generation embedded 
in the arms of a diamond array, the interplay of peri¬ 
odicity along the x-axis and the aperiodic order in the 
transverse directions, produces a highly complex disper¬ 
sion pattern. The flat, dispersionless bands densely fill 
the gaps between the dispersive ones, giving rise to the 
possibility of a quasi-continuous crossover between them 
as the aperiodic components grow in hierarchy. 

The magnetic field piercing the elementary plaquettes 
in each case is shown to control the group velocity of the 
electrons, making them more and more immobile as the 
flux $ —>■ <I>o/2, where $o = hcje is the fundamental flux 
quantum. This is an exemplary case of extreme localiza¬ 
tion induced by the mag netic flux as discussed by Vidal 
and co-workers m- [ij. The magnetic field is shown to 
flip the sign of the effective mass multiple times within 
a single Brillouin zone - a remarkable contrast to the 


ordinary periodic linear lattice. The electron-lattice in¬ 
teraction thus can be fine tuned from outside by selective 
choice of the flux threading a plaquette. 

In what follows, we present our results. In section I we 
work out the basic method of analyzing the FBS in an 
elementary diamond array, and compare the result with 
the existing ones. Section II deals with the Fibonacci- 
diamond chain, and in section III we elaborately discuss 
the fractal-diamond networks. In section IV we explicitly 
discuss the generic diamond loop-array where a magnetic 
field controls the effective mass of the electron making 
its sign periodically flipped. In section V we draw our 
conclusions. 


I. THE HAMILTONIAN AND THE PINNED 
EIGENSTATES 

We refer to Fig. [T] Spinless, non-interacting electrons 
are described using a tight-binding hamiltonian in the 
Wannier basis, viz., 

H = eic\ci -1- tij 

® (b> 

where, Ci is the on-site potential and can assume values 
equal to £3 or £4 depending on the site at the vertex hav¬ 
ing a coordination number 2 = 3 (yellow circles), or in 
the bulk, having a coordination number 2 = 4 (black 
circles). Throughout this paper we shall choose £3 = £4 
just to see the effect of the topology of the lattice alone. 
However, the symbols will be in use to facilitate any dis¬ 
cussion. The nearest neighbor hopping integral tij = t 
along the arm of the diamond, and tij = t along the di¬ 
agonal connecting the vertices with coordination number 
two. The Schrodinger equation, written equivalently in 
the form of the difference equation, 

{E - e^)4’i = ^ Ujil^j (2) 
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allows us to decimate out the “black” vertices of the di¬ 
amond network to map the original array on to an ef¬ 
fective two-arm ladder (Fig. [Ijb)) (with arms marked A 
and B) comprising identical (green colored) atomic sites 
with renormalized on-site potential e = £ 3 -|- 2 t^/(A —£ 4 ). 
The renormalized hopping integral along the arm of the 
ladder now becomes r = t^ j(E — £ 4 ), and the inter-arm 
hopping becomes 7 = t -I- 2t^ /{E — £ 4 ). The decimation 
generates a second neighbor hopping (brown line) inside 
a unit plaquette of the ladder and along the diagonal, 
viz., ^ = t'^/{E - £ 4 ). 

The difference equation for the ladder network may 
now be cast using 2 x 2 matrices, in the form [ 2 ^: 
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i’nA 
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f 4’n+l,A 
\ i^n+l.E 


T ^ 
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'*An—1,A 
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( 3 ) 


It is easy to check that both the ‘potential matrix’ (com¬ 
prising e and 7 ) and the ‘hopping matrix’ (with t and 
commute, and hence can be simultaneously diagonalized 
by a similarity transform. Eq. ([3]) can then be easily de¬ 
coupled, in a new basis defined by The 

matrix M diagonalizes both the ‘potential’ and the ‘hop¬ 
ping’ matrices. The decoupled set of equations are free 
from any ‘cross terms’ and reads, in terms of the original 
on-site potentials and hopping integrals, as: 


E - ea -b t -f 




E — €4 


2t^ 


Pn,A = 


E — €4 


{4'n+l,A + 


(E — €3 + i)(j)n,B — 0 


The first equation represents a periodic array of iden¬ 
tical atomic sites with renormalized on-site potential 
es+t+At^/{E—e 4 ), and nearest neighbor hopping integral 
2t^/{E — € 4 ). The second equation in Eq. (j4|) represents 
(in the new basis) an effective atom, decoupled from its 
neighbors. The potential of this isolated atomic site is 
63 — t. This leads to an eigenfunction with amplitudes 
pinned at the z = 3 vertices as shown in Fig.[T])a) and the 
corresponding eigenvalue is at E = 63 —t. The amplitude 
V'i = 0 at all the z = A vertices for this special energy. 
The non-zero amplitudes are thus trapped in local clus¬ 
ters {z = 3 vertices), as discussed in the introduction. 
The result obtained by Hyrkas et al. @ identifies this 
state as a flat, non-dispersive one (^ . 

To confirm the non-dispersive character of this pinned 
eigenstate, we refer to Fig.[IJc), where an effective linear 
chain of identical atoms is obtained by decimating the ea 
vertices. The resulting on-site potential and the nearest 
neighbor hopping integral for this linear lattice is given 

by, 


Eq — €4 + 


At'^{E - €3 +i) 

{E - £3)2 - <2 


to — 


2f^[E — €3 i) 
{E - £3)2 - <2 


( 5 ) 


ka\ 

{E-e 3 +i) {E - e 4 ){E - € 3 -i) - 8P cos^ { —j =0 

( 6 ) 

In the above equation a is the lattice constant of the 
effective periodic chain. Eq. (I 6 |) clearly indicates a non- 
dispersive, wave vector-independent, flat band at E = 
£3 — t. This is compliant with the result obtained by 
analyzing the decoupled set of equations Eq. (|4]). The 
same set of arcuments brings back the FBS at E = €3 as 
shown in ref [y] for £3 = 0 . 


II. THE FIBONACCI-DIAMOND ARRAY 

We now construct a periodic array of diamond with a 
Fibonacci segment of two different bond lengths L and 
S embedded in it. The linear Fibonacci chain grows ac¬ 
cording to the prescription [ 2 ^ L -A LS and S ^ L. 
This generates a quasiperiodic chain with two differ¬ 
ent hopping integrals, viz., t^ and ts and three kinds 
of vertices a (flanked by two L-bonds), (3 {L on the 
left and S on the right), and 7 (in between an S' — L 
i.A^air). The corresponding on-site potentials are desig¬ 
nated by £ 0 ;, £/3 and £.y respectively. Fig. [2] describes 
periodic array of diamonds with each arm hosting a 
third generation Fibonacci segment LSL. A magnetic 
flux pierces each diamond. The hopping i nteg rals along 
the bonds L and S pick up Peierls’ phases [^, and read 
= tLexp(±z 0 L) and exp(±z 0 s) respec¬ 

tively. Here, 9^s) = 27r$ai(s)/(F„_iaL -I-F’„_ 2 as), En 
being the Fibonacci number in the n-th generation. The 
symbols / and b refer to the hopping in the so called 
forward and backward directions respectively - a conse¬ 
quence of the broken time reversal symmetry. In what 


A (a) 


The linear chain described by Eq. has a dispersion 
relation, 



FIG. 2: (Color online), (a) A perfectly periodic array of iden¬ 
tical diamond blocks when each arm of a diamond hosts a 
third generation Fibonacci sequence of bonds L and S. 

follows, we shall deal with generations ending with the L 

bond only, that is, G' 2 n-i-i, n = 1, 2, . This is only 

for convenience, and does not affect the final result as we 
are interested in the case when n ^ 00 . 

As proposed, we embed a (2n -|- I)-th generation Fi¬ 
bonacci segment in each arm of the diamond. The differ¬ 
ence equations Eq. for the three kinds of atomic sites 
are now of the form, 

{E-CoflPj = tLe^^^ifj+i +tLe~"^^'il;j-i 


{E-ep)ifj = 

{E - e^)ipj = tL( 


'j+i+tLe 


( 7 ) 


where, the index j refers to a site of type a, /3 or 7 
appropriately. 
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The Fibonacci segment trapped in a diamond arm is 
now decimated n times, using the set of Eq. 0 to reduce 
the geometry into a simple diamond with two kinds of 
vertices with renormalized on site potentials (black 
dots) and Cc.n (red dots encircled with green line). The 
effective hopping integral connecting these two sites is 
■> depending on the sense of traversal. The recursion 
relations we exploit during the sequential decimation pro¬ 
cess are given by, 


— ^a,n — l “b 


C/3,n — 


e-y,n — 


Ca,n —1 “b 
^ 7,71 — 1 “b 


^L,n — 


tq = 

S,n 


Cct,n —1 “b 2 
(ff \ 2 -.f 
A„_i 


(e/ 3 ,n-1 + £7,71-1)] 

Ati-I 

{E — f-p^n-l)tL,n-l^L,n-l ^L,n-l^L,n-l 

Ati — I E 

{E — €7,rt-l)^L,77-1^1,77-1 ^5,77-1^5,77-1 

Ati—i E e/ 3 , 77—1 

^5,77-1^1,77-1 [^-E' ~ (e/3,77-1 + £7,77-1)] 


A„- 


77—1 


( 8 ) 


Obviously, ^1(5) „ = ^l(5) n n-th stage of renor¬ 

malization. The quantity A„ is given by A„ = {E — 

^P,n){E — €j^n) — ts,rJ'^S,n- 

As the Fibonacci segment is decimated completely fol¬ 
lowing an n step execution of the recursion relations ([ 5 |) , 
the Fibonacci-diamond geometry is reduced to a simple 
diamond array, as shown in Fig. [T] with ea,n and ec,n 
playing the roles of £3 and €4 in Eq. ([ 5 |) respectively. Ob¬ 
viously, now replaces t and t = 0 . The simple anal¬ 
ogy reveals that, we have localized eigenstates pinned at 
the effective a sites (by virtue of the second equation in 
Eq. (jUl) for all energy eigenvalues obtained by solving 
the equation E — Ca^n = 0 . The corresponding dispersion 
relation, in analogy with Eq. m is given by. 


(E £ 7 : 7 , 77 ) 


(E f^c,n){E £77,77) St^^^cos ( 


ka\ 


= 0 


( 9 ) 

The flat, non-dispersive fc-independent bands are easily 
seen to originate from the solution of the equation E — 

^a,n — 0- 


The solutions of the equation E = ea,n constitute the 
eigenvalue spectrum of a Fibonacci chain as n —>■ 00. 
The spectrum exhibits a global three subband structure. 
Each subband can be scanned over finer scales of the 
wave vector to bring out the inherent self similarity and 
multifractality, the hallmark of the Fibonacci quasicrys¬ 
tals [2^. In the Fibonacci diamond array, we encounter 
precisely this feature, as already evident from Fig. [3l 
Each of the three subbands are populated by the disper¬ 
sive as well as the non-dispersive FBS. The self similarity 


of the bands have been checked by going over to higher 
generations, though we refrain from showing these data 
to save space here. 

In the limit n —>■ 00, that is, when a single diamond arm 
hosts an infinitely large Fibonacci segment, the disper¬ 
sive and non-dispersive bands get more and more densely 
packed, and if one travels along a vertical line at a fixed 
value of the wave vector, a dispersive (non-dispersive) to 
non-dispersive (dispersive) crossover within a single sub¬ 
cluster of states is likely to be observed. Needless to say, 
such a crossover can take place over an arbitrarily small 
interval Afc of the wave vector. In Fig. [3 we have shown 
the bands for a pure transfer model [ 23 |, where the on¬ 
site potentials are set equal to zero, and the quasiperiodic 
order is built in the distribution of the hopping integrals 
only. The central state at E = 0 is there for all values of 
the flux, and as we have checked, for all generations. This 
state belongs to the spectrum of an infinite Fibonacci 
chain, as is well known in the existing literature [2^. 

A magnetic field is found to flatten the dispersive 
bands, decreasing the group velocity, and finally lead¬ 
ing to an extreme localization of the electronic states. 
The spectrum then consists entirely of dispersionless, flat 
bands, grouped in triplicate multifractal families. The 
situation is illustrated in Fig. | 3 ] where, the panels (a) to 
(d) sequentially represent the grouping of flat and disper¬ 
sive bands for a diamond array with each arm hosting a 
5 -th generation Fibonacci chain for <I> = 0 , $ = 4 >o/ 4 , 
4 ) = 24 )o /5 and 4 ) = 4 )o /2 respectively. The gradual 
flattening of the dispersion curves, leading finally to the 
extreme localization in panel (d) is self explanatory. 


III. FLAT BANDS IN THE BERKER-DIAMOND 
ARRAYS 

In this section we discuss two cases where each arm 
of a single diamond hosts hierarchically grown fractals of 
the Berker class. The essential difference with the earlier 
case of a Fibonacci diamond is that, here, with increasing 
hierarchy, the lattice grows in an hierarchical order in the 
transverse direction relative to the arm of a diamond. 
Periodicity is of course, still maintained in the horizontal 
direction. Each unit is thus an approximant of the true, 
infinite fractal. We begin with an open ended Berker- 
diamond array whose growth is illustrated in Fig. SDa). 


A. The open ended Berker geometry 

The linear periodic array of approximants of an open 
ended Berker-diamond array is shown in Fig. IHb). An 
effective two-arm ladder (Fig. HDc)) is generated by dec¬ 
imating out the vertices trapped inside the red dashed 
boxes. The ladder generated in this way develops, by con¬ 
struction, a diagonal (second neighbor) hopping (brown 
line, marked ^) which is equal to the effective nearest 
neighbor hopping r along an arm. The renormalized on- 
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FIG. 3: (Color online) Dispersion curves for a 5th generation Fibonacci-Diamond array for (a) <I> = 0, (b) <I> = $o/4, (c) 
"F = 2>Fo/5, and (d) = "Fo/S. We have chosen Cq, = = 0, tz. = 1 and ts = 2. Energy is measured in units of II- 


site potentials and the hopping integrals of this effective 
ladder geometry are written, using the same symbols as 
in Eq. (| 31 ) as, 


e 


T 


7 


e2,ra 


2Si^nt 


2 

n 



<^2,™ 

2Si^ntn 

^ 2,71 



^ 2,71 


( 10 ) 


where, = {E - e2,n)iE - e3,„) - , and 62,71 = (E - 

e3,n)(<5i,n — tn)- III Ih® above set of equations the z = 2 
and 2 = 3 vertices are assigned the on-site potentials £2 
and £3 respectively, though as before, we shall stick to 
setting £2 = £3. The difference equation (Eq. ([ 3 ])) for 


this new ladder network can again be decoupled, in a 
new basis, to yield the pair of equations. 


E — {^2,71 + 




4 * 71 ,A 


{E ^2,7l)4*7l,B 


-^{4‘71 + 1,A + 4*71 
62 

0 


-i,a) 

( 11 ) 


Naturally, a group of pinned, EBS are obtained as roots 
of the equation E — €2,71 = 0 . In Fig. HKb) we display the 
distribution of amplitudes of the wave function for E = 
£2 = 0 . The amplitudes are non-zero only at the vertices 
of the smallest squares. One square is ‘separated’ from 
its neighboring ones by a vertex at which the amplitude 
is zero. The non-dispersive character of these states can 
be cross-checked, as before, by generating a linear chain 
connecting the green encircled sites in Fig. SKb). The 
dispersion relation becomes, 

{E - £2.„) [5lr7 - < - <5i.n " < cos(fca)] = 0 ( 12 ) 
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(a) 


n= 1 



n = 2 




FIG. 4: (Color online) (a) Construction of an open Berker-diamond chain, (b) Distribution of amplitudes ipi for a flat band at 
E = €2 = 0. We have set t — 1. (c) Conversion of an open Berker-diamond array into a two arm ladder network. 


B. The closed loop Berker geometry 

In a similar manner, a closed loop Berker-diamond ar¬ 
ray, shown in Fig.[5l can be mapped into an effective lad¬ 
der (drawn by blue lines) by decimating out the vertices 
trapped inside the red-dashed boxes. There are two kinds 
of vertices now, viz., with z = 2 and z = 4. The on-site 
potentials bear the symbols €2 and £4 respectively. The 
FBS are extracted by solving the equation E — e 2 ,n = 0. 
In Fig. [5] we demonstrate the distribution of only the flat 
band states against varying magnetic flux for open-ended 
(Fig. EKa)) and closed-loop (Fig. EKb)) Berker-diamond 
arrays. The figure brings out the interesting contrast 
between the two lattices. The FBS for the open ended 
Berker geometry are grossly divided into two segments 
above and below the non-dispersive state at if = 0. The 
subgroups of such states touch each other at if = 0 , 
mimicking a Dirac cone. On the other hand the FBS for 
the closed loop diamond hierarchical array exhibit global 
gaps along the energy axis as well as for changing values 
of the magnetic flux. The edge states in both the cases 
are seen to be most affected by the magnetic field, which 
is likely to lower the persistent current for such energy 
regimes. 

Before ending this section, it is pertinent to comment 
that, in all the cases discussed so far, the amplitudes of 
the FBS are confined within a single unit cell of the lattice 
under consideration. Such states thus belong to the U = 


1 class of FBS as per the nomenclature introduced by 
Danieli et al. Q. However, the rule of constructing the 
distribution of amplitudes, as shown as an example in 
Fig. Hb) can be implemented for any arbitrarily large 
generation of the lattice. This is simple, as any large 
lattice can be renormalized precisely to the one displayed 
in Fig. IlKb). We need to fix the amplitudes of the flat 
band wave function to the values 0, ±1 and ±1/2 on the 
appropriate vertices of the renormalized lattice, exactly 
in the way depicted in Fig. [51 One can then easily work 
in the backward direction to extract the distribution in 
the bare scale of length [ 2 ^ . 

IV. CONTROLLING EFFECTIVE MASS WITH 
MAGNETIC FLUX 

The analyses presented in the previous sections lead to 
an interesting prospect. The effective mass of an electron 
travelling in these kinds of arrays can be controlled by 
tuning the external magnetic field in a non-trivial fashion. 
This result is a generic feature of an array of diamond net¬ 
work, and is thus common to the Fibonacci-diamond and 
Berker-diamond class. The microscopic variations are, of 
course sensitive to the specific aperiodic geometry one 
introduces along the arm of a diamond. We explain the 
basics with the help of the simplest model of a diamond 
array, as depicted in Fig. [I] but now without the diago¬ 
nal hopping t. A uniform magnetic flux $ pierces every 
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FIG. 5: (Color online) (a) Generation and growth of a closed 
loop Berker-diamond array, (b) Distribution of amplitudes ipi 
for a flat band a.t E = 2. We have set t = 1. (c) Conversion 
of a closed loop Berker-diamond array into a two arm ladder 
network. 





FIG. 7: (Color online). Fingerprints of the group velocity of 
the electron in an elementary diamond array. The sections 
between the red lines are filled with a continuous distribution 
of Vg. We set £2 = £4 = 0, t = 1, and display only a few 
contours to help appreciate the discussion in the text. 


(a) 


(b) 






FIG. 6 : (Color online). Distribution of the flat, non-dispersive 
states only for (a) an open ended Berker-diamond array in its 
3rd generation, and (b) a closed loop Berker-diamond array 
in its 3rd generation. We set £2 = £4 = 0, t = 1, and energy 
is measured in units of t. 


plaquette. Such a geometry can easily be mapped on to 
an effectively linear chain, such as depicted in Fig. [IJc). 
The renormalized on-site potentials and nearest neighbor 
hopping integrals on this chain, with £2 = £4 are given 

by, 


t' 


£ + 




E-e 


2t^ cos(7r$/$o) 


(13) 


The dispersive bands are given by, 


E = e±2t 


< 1 -I- cos 



cos(fca) 


(14) 


Before presenting the variation of the effective mass with 
flux, it is essential to remind ourselves that, as an elec¬ 
tron travels around a trapped magnetic flux, the wave 

function picks up a phase, viz., ij} —>■ ■i/io exp j> A.dfj, 
the line integral in the exponent being the flux trapped 
in the closed loop. The magnetic flux here plays a role 
equivalent to the wave vector (2^. One can thus con¬ 
ceive of a k — A> diagram which is equivalent to a typical 
kx — ky diagram for electrons travelling in a two dimen¬ 
sional periodic lattice. The “Brillouin zone” equivalents 
are expected to show up, across which variations of the 
group velocity and the effective mass will take place. 

This is precisely the situation as depicted in Fig. [T] 
Every contour displayed corresponds to a fixed value of 
the group velocity which exhibits a period equal to 2 $o. 
The red lines are the equivalents of the Brillouin zone 
boundaries across which the group velocity flips its sign 
if one moves parallel to the $-axis at a fixed value of fc, or 
vice versa. This implies that, one can, in principle, make 
an electron retard without changing its energy by tuning 
the external flux alone. The group velocity is exactly zero 
along the red lines, showing that the eigenfunctions are 
self-localized around finite clusters of sites, making the 
electronic state a non-dispersive, flat one, as we discussed 
in the beginning. 


E-e 




















However, a more serious issue crops up as we note that 
the same numerical value of the group velocity may lead 
to a positive, or a negative effective mass of the elec¬ 
tron depending on the curvature of the dispersive band 
at those particular values of the group velocity. This is 
not apparent from the fc — $ diagram in Fig. [3 but be¬ 
comes explicit if we look at the expression of the effective 
mass which is given by, 


„ 2 ■^ cos(jS-)cosfca ^ cos2(.|£) sin2(feo) 

^l-|-cos(f|)cosfea 2[l-|-cos(|i) cos(fea)]3/2 

(15) 

In Fig. [ 8 ] we display the variation of the effective mass 
against a changing magnetic field (flux) for different val- 
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FIG. 8: (Color online). Effective mass of an electron travers¬ 
ing a simple diamond array plotted against the magnetic flux 
at (a) ka = 7r/4, (b) ka = 27r/5, (c) ka = 7r/2, and (d) 
ka — 57r/6. We have set e = 0 and t = I and energy is 
measured in unit of t. 


ues of the wave vector k. For 0 < ka < tt/2 , the effective 
mass displays a re-entrant behavior, going from a nega¬ 
tive value below $ < $ 0/2 to a positive value above 
the half flux quantum through a divergence shown pre¬ 
cisely at the half flux quantum. The exception to this 
general character takes place at ka = 7r/2, which is the 
first quarter of the band. The effective mass remains 
negative for all values of the magnetic flux, diverging at 


$ = (p -H l/2)$o, with p = 0, ±1, ±2, . The di¬ 

vergences and crossover in the sign of the effective mass 
are repeated in the interval 1/2 < d>/<i>o < 3/2. The 
variation of m* is flux periodic with a period equal to 
$/$o = 2. 

This discussion points to the obvious implication that 
the electron-lattice interaction can be tuned at will by 
controlling the magnetic flux from outside. The usual 
idea that the electron is made to behave like a positively 
charged particle in a certain part of the band where its 
effective mass turns negative yields, in this case, a wider 
variety of phenomenon where this can be achieved in a 
re-entrant fashion throughout the band, by tuning an 
external perturbation, that is, the magnetic field. To 
the best of our knowledge, this simple aspect of such 
problems was not highlighted before. 


V. CONCLUSIONS 


We have shown how to extract the eigenvalues corre¬ 
sponding to the localized, non-dispersive, degenerate flat 
band eigenstates of an infinite, periodic diamond loop 
array, where aperiodic structures of increasing hierarchy 
grow on each arm of a basic diamond network. A real 
space renormalization group scheme is exploited to un¬ 
ravel a countable infinity of such states. The grouping of 
the dispersive and non-dispersive bands in each case re¬ 
sembles the actual band structure of the quasiperiodic or 
fractal lattices considered here, as their generation tends 
to the respective thermodynamic limits. An external 
magnetic field is shown to be able to control the cur¬ 
vature of the dispersive bands, and hence the sharpness 
of localization. The group velocity and the effective mass 
of the electron are thus shown to exhibit a re-entrant be¬ 
havior inside a single Brillouin zone, a fact that turns 
out to be a consequence of the quasi-one dimensionality 
of the loop structure. The scheme is easily extendable to 
photonic, phononic of magnonic excitations, and the flat 
band eigenstates for an aperiodically grown superlattice 
can thus be worked out. This may throw new challenges 
to experimentalists to engineer the formation and posi¬ 
tioning of the non-dispersive energy bands in such artifi¬ 
cial lattices. A possible application in device technology 
may thus be on the cards. 
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